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ABSTRACT 

We present our successful program using Chandra for identifying the X-ray afterglow with 
sub-arcsecond accuracy for the short GRB 111117A discovered by Swift and Fermi. Thanks to 
our rapid target of opportunity request, Chandra clearly detected the X-ray afterglow, whereas 
no optical afterglow was found in deep optical observations. Instead, we clearly detect the host 
galaxy in optical and also in near-infrared bands. We found that the best photometric redshift 
fitof the host is z == 1.3lio;^f (90% confidence), making it one of the highest redshift short GRBs. 
Furthermore, we see an offset of 1.0±0.2 arcseconds, which corresponds to 8.4± 1.7 kpc assuming 
z=l .31, between the host and the afterglow position. We discuss the importance of using Chandra 
for obtaining sub-arcsecond localization of the afterglow in X-rays for short GRBs to study GRB 
environments in great detail. 

Subject headings: gamma rays: bursts 
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1. Introduction 

Gairima-ray bursts (GRBs) are traditionally di- 
vided in two classes based on their duration and 
spectral hardness: the long duration/soft spec- 
trum GRBs, and the short duration/hard spec- 
trum GRBs (Kouveliotou et al. 1993). The two 
classes of bursts further differ in their spectral 
lags, the measurement of the delay in the arrival 
time of the low-energy photons with respect to 
the higher energy ones: long bursts tend to have 
large positive lags, while short bursts exhibit neg- 
ligible or negative lags (Norris & Bonnell 2006). 
The long standing paradigm is that these two phe- 
nomenological classes of GRBs originate from dif- 
ferent progenitor systems. A preponderance of 
evidence now links long GRBs with the death of 
massive stars (Woosley & Bloom 2006, and refer- 
ences therein), yet the origin of short GRBs re- 
mains largely unknown. The common notion that 
short bursts originate from coalescing compact bi- 
naries, either neutron star-neutron star (NS-NS) 
or neutron star-black hole (NS-BH) mergers (e.g., 
Eichler et al. 1989; Paczynski 1991; Narayan et 
al. 1992; Rosswog 2005; Rezzolla et al. 2011), 
makes them the most promising tool to aid in 
the direct detection of gravitational waves (GWs) 
by forthcoming facilities such as Advanced-LIGO, 
Advanced- VIRGO or LCGT (KAGRA) (e.g., Nis- 
sanke et al. 2010). It is therefore of primary im- 
portance to convincingly corroborate the mergers 
scenario with a robust observational basis. 

Significant progress in understanding the ori- 
gin of short GRBs has been achieved only re- 
cently. This advance was enabled by the detec- 
tion of their afterglows in 2005 thanks to the rapid 
position notice and response by HETE-2 (Ricker 
et al. 2003) and Swift (Gehrels et al. 2004). 
The first X-ray afterglow of a short GRB was dis- 
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covered in GRB 050509B by Swift (Gehrels et ad. 
2005; Bloom et al. 2006), allowimg us to localize a 
short GRB within a ~5" accuracy, to identify its 
host galaxy and therefore determine its redshift. 
Just two months after this discovery, the first opti- 
cal afterglow of a short GRB was discovered from 
GRB 050709 detected by HETE-2 (Villasenor et 
al. 2005; Fox et al. 2005; Hjorth et al. 2005). 
The. optical position, which is characterized by a 
sub-arcsec accuracy, allowed us to unambigously 
identify the GRB host, but also to precisely pin- 
point the GRB location within its galaxy, in a re- 
gion of low-star formation. The very first local- 
izations of short GRBs immediately provided us 
with fundamental clues about their nature. They 
demonstrated that short GRBs are cosmological 
events with an energy scale of 10 49 -10 51 erg, that 
they occur in a different environment from long 
GRBs, and are not associated with bright Type Ic 
supernovae (Bloom et al. 2006; Prochaska et al. 
2006; Covino et al. 2006). 

Since 2005 the sample of well-localized short 
GRBs has significantly grown, allowing for a 
deeper insight into the nature of their progeni- 
tors. The observed redshift distribution, ranging 
0.11 < z < 1, hints at a progenitor system with a 
broad range of lifetimes (Berger et al. 2007) . This 
is also consistent with the heterogeneous sample of 
galaxies harboring short bursts (e.g. Malesani et 
al. 2007; D’Avanzo et al. 2009; Leibler & Berger 
2010). According to Berger (2011), ~50% of 
short GRBs occur in star forming galaxies, 10% of 
short GRBs are associated with elliptical galaxies, 
while the remaining 40% cannot be classified due 
to various observational constraints. Although a 
larger fraction of short GRB hosts are star forming 
galaxies, they differ from star forming host galax- 
ies of long GRBs, which have lower masses and 
younger population ages (Gorosabel et al. 2006; 
Leibler & Berger 2010). Another critical test of 
the progenitor models is the observed offset dis- 
tribution of short bursts (TVoja et al. 2008; Fong 
et al. 2010; Church et al. 2011). The median 
physical projected offset between the host center 
and the short GRB position is ~5 kpc (Fong et al. 
2010), which is about five times larger than that 
of long GRBs (Bloom et al. 2002), and shows a 
broader dispersion. This is in agreement with the 
merger scenario, as in several models NS-NS/NS- 
BH systems are expected to receive significant kick 
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velocities at birth (Bloom et al. 1999; Fryer et al. 
1999; Belczynski et al. 2006), or to dynamically 
form in globular clusters in the outskirts of their 
galaxies (Grindlay et al. 2006). 

Despite the major progress of the last few years, 
the study of short GRBs and their progenitors 
has still been suffering from their less secure after- 
glow positions and redshifts. Unlike long GRBs, 
none of the redshifts of short GRBs has been di- 
rectly measured through afterglow spectroscopy, 
and only in the case of GRB 060121, a photomet- 
ric redshift was derived from the afterglow spec- 
tral energy distribution (de Ugarte Postigo et al. 
2006). This is because the optical afterglows are 
significantly fainter than those of long GRBs (Ny- 
sewander et al. 2009; Kann et al. 2011). Almost 
all of the redshifts of short GRBs are instead mea- 
sured from spectroscopic observations of the “as- 
sociated” host galaxy. The likelihood of a spurious 
association is small when a sub-arcsec position is 
available (Cobb et al. 2006). However, if an after- 
glow is only detected by the Swift X-Ray Telescope 
(XRT; Burrows et al. 2005), a chance alignment 
becomes more significant due to the larger uncer- 
tainty in the localization (2-5"). Unfortunately, 
the latter scenario represents the majority of cases 
(~65%). 

A further bias is introduced by the fact that 
sub-arcsecond positions are mainly derived from 
optical afterglow detections. In the standard fire- 
ball model (Sari et al. 1998), the optical bright- 
ness depends sensitively on the density of the cir- 
cumburst environment. This effect disfavors the 
accurate localization of short GRBs occurring in 
the lower-density galaxy halo or even outside their 
own galaxy, in the intergalaCtic medium. Such 
populations of large-offest short GRBs has already 
been hinted at Bloom et al. (2007) and TYoja 
et al. (2008). However being localized mainly 
by XRT, their association with the putative host 
galaxy remain uncertain. Increasing the sample of 
large-offset short bursts with sub-arcsecond local- 
ization is crucial to discriminate whether their pro- 
genitors were ejected from their birth site, favoring 
models which predict NS binaries with large kick 
velocities and ~Gyr lifetimes, or they were formed 
from dynamical interactions in globular clusters 
(Salvaterra et al. 2010). 

In this context, rapid Chandra observations of 
short GRB afterglows represent the critical ob- 


servational gateway to overcome the current ob- 
servational limits. Since 80% of the Swift short 
GRBs with prompt slewing are detected in X-rays, 
and only 40% of them are detected in the optical 
band, X-ray observations have a higher probabil- 
ity of successfully detecting the afterglow of short 
GRBs. The superb angular resolution of Chandra 
allows for a sub-arcsecond localization compara- 
ble to the optical ones, thus enabling the secure 
host identification and the precise measurement 
of the GRB projected offset. Furthermore, be- 
cause the X-ray afterglow is less subject to density 
effects, Chandra localizations allow us to build a 
sample of well-localized short GRBs with limited 
bias. This is the key to distinguish among possi- 
ble different short GRB populations (Sakamoto & 
Gehrels 2009), which couid arise from a different 
progenitor and/or environment. 

In this paper, we report the first results of our 
Chandra program which led to the accurate local- 
ization of GRB 111117A detected by Swift and 
Fermi. Our results were leveraged with an in- 
tense ground-based follow-up campaign. No opti- 
cal/infrared counterpart was found, therefore our 
Chandra localization uniquely provides the most 
accurate sub-arcsecond position. The paper is or- 
ganized as follows: we report the prompt emission 
properties in §4, the X-ray afterglow properties in 
§5.1, the deep optical afterglow limits in §5.2, and 
the host galaxy properties in §6. We discuss and 
summarize our result in §7. The quoted errors 
are at the 90% confidence level for prompt emis- 
sion and X-ray afterglow data, and at the 68% 
confidence level for optical and near infrared data 
unless stated otherwise. The reported optical and 
near infrared magnitudes are in the Vega system 
unless stated otherwise. Throughout the paper, 
we use the cosmological parameters, fl rn = 0.27, 
Ga = 0.73 and Ha = 71 km s -1 Mpc -1 . 

2. GRB 111117A 

On 2011 November 17 at 12:13:41.921 UT, the 
Burst Alert Telescope (BAT; Barthelmy et al. 
2005a) triggered and localized short GRB 111117A 
(Mangano et al. 2011). The Fermi Gamma-ray 
Burst Monitor (GBM; Meegan et al. 2009) also 
detected the burst (Foley et al. 2011). The 
BAT location derived from the ground analysis 
was (R.A., Dec.) (J2000) = (00 h 50 m 49.4 s , +23° 
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00 ; 36") with a 90% error radius of 1.8'. The Swift 
XRT started its observation 76.8 s after the trig- 
ger. A fading X-ray source was found at the lo- 
cation of (R.A., Dec.) (J2000) = (00 h 50 m 46.22 s 
, +23° 00' 39.2") with a 90% error radius of 2.1" 
(Melandri et al. 2011a). The Swift UV-Opticai 
Telescope (UVOT; Roming et al. 2005) began 
the observations of the field 137 s after the trig- 
ger, and no optical afterglow was detected (Oates 
et al 2011). 

The earliest ground observations of the field 
were performed by the Gao-Mei-Gu telescope 
(GMG) at 1.96 hr after the BAT trigger. No after- 
glow was detected with the exposure time of 600 s 
'within the XRT error circle by GMG in the R band 
(Zhao et al. 2011). The Nordic Optical Telescope 
(NOT) observed the field at 8.9 hr after the burst, 
and found an optical source inside the XRT error 
circle (Andersen et al. 2011). The optical source 
wa3 confirmed with a possible extended morphol- 
ogy by the Magellan/Baade telescope (Pong et 
al. 2011), the Gemini-South telescope (Cucchiara 
et al. 2011), the GROND telescope (Schmidl et 
al. 2011), and the Telescopio Nazionale Galileo 
(TNG; Melandri et al. 2011b). The Gran Tele- 
scopio CANARIAS (GTC), the Subaru telescope, 
the United Kingdom Infrared Telescope (UKIRT) 
and the CanadarFrance-Hawai Telescope (CFHT) 
also collected images of the field. 

Based on no clear detection of an optical af- 
terglow of the short GRB 111117A, we triggered 
our Chandra Target of Opportunity (ToO) obser- 
vation 6 hours after the trigger (Sakamoto et al. 
2011b), and the observation started 3 days later. 
The X-ray afterglow was clearly detected in 20 ks, 
obtaining a sub-arcsecond position of the after- 
glow in X-rays (Sakamoto et al. 2011c). 

3. Data analysis 

HEAsoft version 6.11 and the Swift CALDB 
(version 20090130) were used for the Swift BAT 
data analysis. The XRT data products were ob- 
tained from the automated results available from 
the UK Swift Science Data Center (Evans et al. 
2007, 2009). CIAO 4.3 and CALDB 4.4.6 were 
used for the Chandra data analysis. The Fermi 
Gamma-ray Burst Monitor (GBM) data were pre- 


pared using the RMFIT software package, 1 with 
data from three Sodium Iodide (Nal) scintillation 
detectors (detector ID 6, 7 and 9) and two Bismuth 
Germanate (BGO) scintillation detectors (detec- 
tor ID 0 and 1). 

A standard data reduction of optical and near 
infrared images was performed using the IRAF 2 
software package. The SExtractor 3 (Bertin & 
Arnouts 1996), the SkyCat GAIA 4 and IRAF 
were used to extract sources and perform the pho- 
tometry. To accomplish consistent photometry for 
images collected by various telescopes, we selected 
10 common stars in the field and perform a rela- 
tive photometry (Table 1). When some of the stars 
were saturated (especially for a large aperture tele- 
scope such as GTC), a subset of these 10 reference 
stars were used. The USNO B-l R2 magnitude or 
the SDSS magnitudes were used as the reference 
magnitude for the stars. For the near infrared im- 
ages of UKIRT and CFHT, we use the reference 
stars in the 2MASS catalog. The galactic extinc- 
tion has been corrected using E(B-V) = 0.03 mag 
toward the direction to this burst (Schlegel et al. 
1998). The log of optical and near infrared obser- 
vations presented in this paper is summarized in 
Table 2. 

4. Prompt Emission 

The light curve of the prompt emission is com- 
posed of two episodes: the first episode shows 
multiple overlapping pulses with a total duration 
of 0.3 s, and the second episode is composed of 
two pulses with a duration of 0.1 s (Figure 1). 
The duration is Tgo = 464 ± 54 ms (1 a error; 
15-350 keV) measured by the BAT background- 
subtracted light curve using the mask modulation 
(e.g., mask-weighted light curve). This Tgo du- 
ration is significantly shorter than 2 s, which is 
the standard classification of short GRBs in the 
BATSE GRBs (Kouveliotou et al. 1993). Fur- 
thermore, this duration is shorter than 0.7 s, which 
is claimed to be the dividing line between long 
and short GRBs for the Swift sample (Bromberg 
et al. 2012). The hard-to-soft spectral evolution 
is seen in both the first and the second episodes 


1 http://fermi.gs fc.nasa.gov/ssc/data/analysis/user/ 
2 bttp: / / iraf.noao.edu / 

3 http: //www.astromatic.net /software/sextr actor 
4 http: //astro.dur. ac.uk/" pdraper/gaia/ gaia.html 
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of GRB111117A (see the hardness ratio plot at 
the bottom panel of Figure 1). There is no indi- 
cation of extended emission (Norris et al. 2011) 
down to a flux level of ~2 x 10 -10 erg cm“ 2 s -1 , 
assuming a power-law spectrum with a photon in- 
dex of —2 in the 14-200 keV band by examining 
the BAT sky image from 60 s (after the spacecraft 
slew settled) to 950 s after the BAT trigger time 
(hereafter to, bat)' The spectral lag between the 
100-350 keV and the 25-50 keV band is 0.6 ± 2.4 
ms, which is consistent with zero, using the BAT 
raw light curves (non mask-weighted light curves)' 
by subtracting the background by a constant rate 
around the burst. In the fluence ratio versus Tgo 
plane, GRB 111117A is located in the same region 
where most of the BAT short GRBs are populated 
(Figure 2). 

The time-integrated spectral properties are in- 
vestigated by performing a joint spectral analysis 
with BAT and GBM data. The spectrum is ex- 
tracted from to, bat + 0.024 s to to, bat + 0.520 
s using batbinevt for the BAT data and using 
the RMFIT software package for the GBM data 
in the same time interval. The BAT energy re- 
sponse file is generated by batdrmgen. The energy 
response files located at the HEASARC Fermi 
archive for trigger bnlllll7510 are used for the 
GBM. We use the xspec spectral fitting package 
to do the joint fit. The energy ranges of 15-150 
keY, 8-900 keV and 0.2-45 MeV are used for the 
BAT, the GBM-Nal and the GBM-BGO instru- 
ments, respectively. A constant factor is multi- 
plied by the model to take into account the cal- 
ibration uncertainty among the different instru- 
ments. The best fit spectral parameters are sum- 
marized in Table 3. We find that a power- law 
multiplied by an exponential cutoff (CPL) 5 pro- 
vides the best representative model of the data. 
The best fit parameters in this model are the 
power-law photon index o CPL = — 0.28to|g and 
Ep ea k = 420i^ keV (xVd.o.f. = 627/661). The 
90% confidence interval of a constant factor of the 
GBM detectors is between 0.50 and 0.78 which is 
an acceptable range taking into account a current 
spectral calibration uncertainty between the BAT 
and the GBM. A simple power-law model shows 

5 /(E) oc £“ CPL exp ^ E ^ , where a CPL is the 

power-law photon index and E prJl k is the peak energy in 

the v F„ spectrum. 


a significantly worse fit to the data (y 2 /d.o.f. = 
729/662). There is no significant improvement of 
X 2 using a Band function (Band et al. 1993) fit 
(x 2 /d.o.f. = 627/660) over a CPL fit. A pref- 
erential fit to a CPL model and also a system- 
atically harder photon index than that of a long 
GRB seen in GRB 111117A are general character- 
istics of a time-integrated spectrum of short GRBs 
(e.g., Ghirlanda et al. 2009; Ohno et al. 2008). 
The fluence in the 8-1000 keV band calculated us- 
ing the best fit time-integrated spectral parame- 
ters baaed on a CPL fit above is 7.3t 2 '® * 10” 7 erg 
cm -2 . Due to a poor statistics to extract a spec- 
trum in a very short time window, a peak flux 
is calculated by scaling the BAT mask-weighted 
count-rate into a flux by folding the BAT energy 
response and assuming the best fit time- integrated 
spectral parameters in a CPL model. The peak en- 
ergy flux at the 8-1000 keV band in 50 ms window 
starting from to, bat + 0.450 s is (3.8 ± 1.2) x 10~ 6 
erg cm -2 s -1 . The time-resolved spectroscopy is 
difficult to perform due to the limited statistics in 
the data. 

We search for pre-burst emission analyzing 
the BAT survey data (Detector Plane Histogram; 
DPH). Approximately 4.5 hours before the burst 
trigger, GRB111117A was in the field of view 
of BAT (26.1° from the boresight direction) for 
~1 ks during the observation of the blazar PKS 
0235+16 (observation ID 00030880085). We use 
bat survey script to process the DPH data. The 
extracted rates at the location of GRB 111117A 
are corrected to the on-axis rate by applying an 
off-axis correction based on observations of the 
Crab. We find no significant emission during this 
observation at the location of GRB 111117A. The 
3 sigma upper limit assuming a power-law spec- 
trum of a photon index of —2 is 1.4 x 10 -9 erg 
cm -2 s -1 at the 14-200 keV band in 300 s expo- 
sure. 

5. Afterglow 
5.1. X-rays 

The Swift XRT X-ray afterglow light curve can 
be fit to a simple power-law decay (Figure 4). The 
spectrum collected in the Photon Counting (PC) 
mode is well described by an absorbed power- law 
model. The best fit spectral parameters are a 
photon index of -2.19to;3g an< * && excess N a of 
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* 10 21 cm 2 assuming the galactic Nh at 
the burst direction of 3.7 x 10 20 cm -2 . Both the 
measured photon index and Njj of GRB 111117A 
are consistent with other Swift short GRBs. 

The Chandra observation started at 12:39:25 
UT, and ended at 18:39:10 UT on 2011 November 
20 with a total exposure of 19.8 ks. The ACIS 
instrument had five CCD chips (S3, S4, S5, 12 
and 13) enabled, with the S3 chip as the aiming 
point for the source. The data were collected in 
the FAINT mode. The X-ray afterglow is clearly 
detected in the Chandra observation with 3.9 a 
significance by wavdetect (source net counts of 8) 
within the XRT error circle. To refine the astrom- 
etry of the Chandra data, we apply the same anal- 
ysis method described in Feng & Kaaret (2008). 
We extract the Chandra image (0.35 - 8 keV) 
that overlaps with the GTC image (4.4' x 8.7'). 
The astrometry of the GTC image is calibrated 
against the SDSS catalog, and its standard devi- 
ation is ~0.3". We rim wavdetect with options 
of scales— ’1.0 2.0 4.0 8.0 16.0” and sigthresh = 
4 x 10~ 6 to the extracted Chandra image. There 
are four sources which have a good match between 
the images. Then, we use geomap task in the IRAF 
IMMATCH package to find the best coordinate 
transformation between the Chandra and the GTC 
image by fitting those four sources. Finally, we 
apply geoxytrans task (IRAF IMMATCH pack- 
age) for the originally detected Chandra position 
using the coordinate transformation calculated by 
geemap to find the astrometrically corrected Chan- 
dra afterglow position. The refined afterglow po- 
sition is shifted £R.A. = —0.221" and <5Dec. = 
—0 020" from the position originally derived by 
wavdetect. The best Chandra X-ray afterglow po- 
sition is (R.A., Dec.) (J2000) = (00 h 50 m 46.264 a , 
+23° 00' 39.98") with 1 a statistical uncertainty 
of 0.09" in right ascension and 0.16" in declination 
coordinate. The Chandra position is well within 
the XRT 90% error circle (see Figure 3). 

The combined Swift XRT and Chandra X-ray 
afterglow light curve is well fit by a simple power- 
law with the index of — 1.25ionl' As shown in 
Figure 4, the X-ray afterglow of GRB 111 117A be- 
longs to a dim population of the Swift short GRBs. 
When applying the classification scheme of X-ray 
afterglow of short GRBs suggested by Sakamoto 
& Gehrels (2009), the afterglow is classified as a 
’’short-lived” X-ray afterglow since the X-ray flux 


at 10 4 s after the trigger is < 10 13 erg cm 2 s -1 . 

5.2. Optical 

We investigate the possible optical afterglow 
emission by using the image subtraction technique 
between the early and the late time epoch obser- 
vations by TNG and GTC. We use ISIS software 
package (Alard & Lupton 1998) to perform the 
image subtraction. The early and the late epoch 
observations of TNG and GTC were obtained at 
*o,bat + 7.23 hours and to, bat + 7.89 hours, 
and to, bat + 11.4 days and to, bat + 14.4 days, 
respectively. We find no significant emission at 
the Chandra X-ray afterglow location in the sub- 
tracted images in both the TNG and the GTC 
observations (Figure 5), with 3 sigma upper lim- 
its of R > 24.1 mag for TNG and r > 24.9 (AB) 
mag for GTC. The TNG limiting magnitudes of 
the first and second epochs are R > 24.7 mag and 
R > 25.4 mag, respectively. For GTC, the limit- 
ing magnitude of the first and second epoch are 
r > 25.8 (AB) mag and r > 26.1 (AB) mag, re- 
spectively. Those limits are some of the deepest 
optical limits on short GRBs ever obtained (see 
upper panel of Figure 6). 

6. Host Galaxy 

The host galaxy of GRB 111117A has been de- 
tected from the near infrared to the optical bands. 
There is only one near-infrared and optical source 
located near the Chandra X-ray afterglow position. 
Although the weak nature of the source makes it 
difficult to investigate whether the source is ex- 
tended or not, the optical flux of the source is 
constant between 7 hours and 14 days after the 
burst at a level of ~1.1 /aJy (bottom panel of Fig- 
ure 6). Therefore, we conclude that the source 
detected in K, J, z, i, r, g , and R bands is the 
host galaxy of GRB 111117A (Figure 7). 

To estimate the redshift of the host galaxy, we 
perform a spectral energy density (SED) fit with 
the stellar population model of Maraston (2005). 
We use the Single Stellar Populations (SSP) mod- 
els with a Salpeter initial mass function (Salpeter 
1955), solar metallicity which ranges from 0.005 
Z 0 to 3.5 Z G (0.005, 0.02, 0.5, 1.0, 2.0 and 3.5 
Zq), and a red or blue horizontal branch mor- 
phology. A total number of 269 SED templates 
ranging in stellar age from 10 Myr to 15 Gyr was 
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applied. The Bayesian Photometric Redshift soft- 
ware (BPZ; Benitez 2000) is used to fit the data 
in g, r, and i bands (GTC), z band (GMG), J 
band (CFHT) and K band (UKIRT) with those 
SED templates. We found that the best fit SED 
template corresponds to a solar metalicity, a red 
horizontal branch morphology and the luminosity 
weighted mean stellar age of 0.1 Gyr with a red- 
shift of 1.36*5 .is- Our best fit SED template of 
SSP model with a solar metallicity and a red hori- 
zontal branch matches well with other short GRB 
hosts studied by Leibler k Berger (2010). As seen 
in Figure 9 (top), there is a less significant low red- 
shift solution (z < 0.25) in the estimated redshift. 
We found that this low redshift solution is coming 
from the template with the young stellar ages of 
~10 Myr. As we will discuss in §7, it is unlikely 
that the host galaxy has the stellar ages of ~10 
Myr. Therefore, to constrain the redshift better, 
we focus on 41 SED templates of a solar metallic- 
ity and a red horizontal branch morphology with 
stellar age from 20 Myr to 15 Gyr. Figure 8 shows 
the best fit SED and the photometric data. The 
GTC spectrum with the exposure time of 4 x 1800 
s is also shown in the figure. The signal-to-noise is 
low, and there are no clear absorption or emission 
line features in the spectrum. The continuum part 
of the spectrum is consistent with the best fit SED 
template. Bottom panel of Figure 9 shows the 
posterior probability distribution of the estimated 
redshift for this SED fit. No low redshift solu- 
tion is evident in the probability distribution. We 
find the best estimated redshift to be 1.31 (90% 
confidence interval 1.08 < z p h. < 1-77). The likeli- 
hood that the redshift is correct is 80% (reduced 
X 2 of the fit is 0.65). The best fit SED template is 
the case with the luminosity weighted mean stellar 
age of 0.1 Gyr and a mass of ~ 1 x 10 9 M 0 . To 
investigate the likelihood of the host being a star- 
forming galaxy, we also examine the SED tem- 
plates with an exponentially decaying star forma- 
tion rate (Maraston 2010), with an e-folding time 
of 0.1, 1 and 10 Gyr (stellar age ranges from 10 
Myr to 15 Gyr). Although our J band data point 
shows a relatively poor agreement with the best fit 
template with an e-folding time of 0.1 Gyr and the 
luminosity weighted mean stellar age of 0.3 Gyr, 
the fit is still acceptable (reduced x 2 = 0.91). The 
best fit redshift in this case is 1.18t 9 ;fi. The fit 
becomes worse if the e-folding time gets larger. 


Therefore, our current data also supports the host 
galaxy to be ~0.1 Gyr post-star burst galaxy. We 
also fitted the optical-NIR SED using MAGPHYS 
package (de Cunha et al. 2008) to check the va- 
lidity of the fitting result. The MAGPHYS fit in- 
cludes an extinction parameter as a part of the 
fit, but performs the fitting at a given redshift. 
An exponentially decaying star formation rate is 
assumed. The redshifts were increased by a step of 
0.05 from 0.8 to 1.5. We confirm that the returned 
X 2 is the least in the range of the z P h a t from the 
BPZ, and a moderate extinction of Ay = 0.2 — 0.5 
mag is found. The best fit solutions give the ex- 
ponential time scale of about 1.5 Gyr, with the 
luminosity weighted mean stellar age (in r-band) 
of a few hunder Myr, the stellar mass of a few 
times 10 9 M©. These output values are consistent 
with the solutions we got already through during 
the photometric redshift derivation. In summary, 
based on various SED template fits, we can con- 
clude the followings about the host galaxy. First, 
the redshift of the host is ~1.3 regardless of the 
SED model. Second, the host is either a star- 
forming galaxy of the luminosity weighted mean 
stellar age of 0.1 Gyr and a mass of ~ 1 x 10 s 
M 0 or a post-starburst galaxy. Further deep J 
or Y band data are crucial to pin-down the host 
properties. 

A significant offset between the center of the 
host galaxy and the X-ray afterglow has been 
found for GRB111117A. The center position of 
the host galaxy has been examined by running 
SExtractor to the 2nd epoch of the GTC r im- 
age which is our highest quality optical image. 
The best location of the host center is (R.A., 
Dec.) (J2000) = (00 h 50 ra 46.258 s , +23° 00' 
40.97"). The position changes less than half a 
pixel (< 0-13") by changing the detection thresh- 
old of SExtractor. Therefore, the projected offset 
between the center of the host galaxy and the X- 
ray afterglow is 1.0" (dR.A. = 0.083" and <5Dec. 
= —0.990"; see Figure 7). Taking into account the 
statistical error in the X-ray afterglow position of 
0.18" and the statistical error of the host center 
location of 0.13", we estimate the offset with its 
error to be 1.0±0.2 arcseconds, which correspond- 
ing to a distance of 8.4 ± 1.7 kpc at a redshift of z 
= 1.31. This offset is consistent with the median 
offset of 5 kpc for short GRBs (Fong et al. 2010). 

Throughout this section, the association be- 
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tween the GRB and the galaxy close to the Chan- 
dra position as the host has been assumed. It 
is however worth discussing the chance associa^ 
tion probabilities. Using the formula provided by 
Bloom et al (2002), the probability of finding an 
unrelated galaxy of the R magnitude of ~23 .3 with 
the distance of 1.0" is 0.8%. Although this prob- 
ability is comfortably small, it is non-negligible. 
Therefore, we would like to note the small pos- 
sibility of a mis-identification of the host galaxy 
based on our currently available dataset. 


7. Discussion 


Our photometric redshift of 1 - 3115 x 23 (90% con- 
fidence) of the host galaxv is a reasonable esti- 
mate for the following reasons. First, by plotting 
the observed r AB magnitude {tab) of the host 
galaxies of short GRBs as a function of redshifts 
(Berger 2009), we find that the relatively low 
magnitude of the host galaxy of GRB111117A, 
ta 3 = 24.20±0.07, is located at the redshift range 
of > 0.5 (Figure 10). Second, we found that the 
less significant low redshift solution (z < 0.25) in 
the photometric redshift estimation (Figure 9) is 
coming from the templates with the young stel- 
lar ages of ~10 Myr. At z=0.25, if the galaxy is 
star forming, there would be a chance of seeing 
emissions of [Oil], H-/3 and [OHI] in the optical 
spectrum. However, we see none of those lines 
in 'he GTC spectrum. Furthermore, ~10 Myr is 
in general too young for a galaxy, and also is in 
agreement with the lifespan of a 40 Mg star. It is 
very unlikely that NS-NS merger could happen in 
a ~10 Myr-old galaxy. Therefore, the low redshift 
solution for the photometric redshift is unlikely 
to be the case of GRB 111117A. Therefore, here- 
after, we will discuss the rest-frame properties of 
GRB 111117A using our best photometric redshift 
of 1.31. 


Assuming the redshift of 1.31, the isotropic 
equivalent 7 -ray energy (£ 7 ,iso) which is inte- 
grated from 1 keV to 10 MeV at the rest-frame 
is 3 . 4^5 x 10 51 erg. The peak energy at the rest- 
frame (Fp^) is 945l|jQ keV. The 90% errors in 
£7, iso and E^ ?k are taking into account not only a 
statistical error but also an uncertainty in the es- 
timated redshift. As shown in Figure 11, F 7i ; m of 
GRB 111117A is located at the high end of E lt \ so 
distribution of short GRBs and at the low end 


of £ 7 , iso distribution of long GRBs. Relatively 
low £ 7 , ^ and high E^ k compared to those of 
long GRBs make GRB 111117A inconsistent with 
the £^,-£ 7,180 (Amati) relation (Amati et al. 
2002). This characteristic is consistent with being 
a short GRB because most of the short GRBs are 
well known outliers of the Amati relation (Amati 
2006; Nava et al. 2012). Recently, Margutti et al. 
(2012) and Bernardini et al. (2012) found a uni- 
versal scaling among £ 7,1301 £peak and the X-ray 
energy emitted in the 0.3-10 keV band (£x,iao)- 
This three-parameter relation, the £x,iso-£ 7 ,iso- 
£pg C ak relation, is claimed to hold not only for 
long and short GRBs, but also for low luminos- 
ity GRBs. The best fit function is log(£x,iso) = 
1.061og(£7, i8O ) - 0.741og(£^ k ) - 2.36, with the 
intrinsic scatter of 0.31 (1 a). It is interesting to 
see whether GRB111117A is consistent with the 
£x , iso- £ 7 .iso- £peak relation or not. By integrat- 
ing the best fit power-law function of the X-ray 
light curve data of XRT and Chandra from the 
start time (to, bat + 90 s) to the end time (to, bat 
+ 282 ks), the X-ray fluence in the 0.3-10 keV 
band is estimated to be 2.3 x 10 -8 erg cm -2 . Us- 
ing the redshift of 1.31, £x,i SO of GRB111117A 
is calculated to be 1.0 x 10 so erg. Therefore, 
GRB111117A is located well within 2-<r region 
from the best fit £x,:so-£ 7 ,iscr £p r e c ak relation. Fur- 
thermore, GRB 111117A is populated at the same 
region where seven short GRBs are located in 
the £x,iso-£ 7 ,iBo-£p r e c a j t relation (see Figure 20 of 
Margutti et al. 2012). 

The upper limit on the afterglow emission in 
optical and the weak X-ray afterglow emission 
are consistent with a standard afterglow model if 
a cooling frequency (u c ) is located above both 
the optical and the X-ray band with a rela- 
tively shallow power-law index (p) for the elec- 
tron distribution (p < 2.1), The optical-to-X- 
ray spectral index (Jakobsson et al. 2004), /fox 
(= log{/^(£)/ fv (3 keV)}/ log(iy 3 keV / im)) , is esti- 
mated to be < 0.78 using the same definition on 
the X-ray flux density measured at 11 hours after 
the burst at 3 keV, and the optical afterglow limit 
based on the GTC r band. This upper limit of 
/fox is within the allowed range of the standard 
afterglow model between 0.5 to 1.25. Following 
the discussion on Nysewander et al. (2009), 
we calculate the flux ratio between the optical 
( Fr ) and the X-ray (Ex) to constrain the stan- 
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dard afterglow emission model. The flux ratio of 
GRB111117A is estimated to be Fr/Fx < 590 
using the X-ray flux at 11 hours after the burst at 
5 keV, which is the same definition of Nysewander 
et al. (2009), and the upper limit in r band of 
GTC. If a cooling frequency is above the optical 
and the X-ray range (i/ c > vx > i'r), the ratio is 
expressed as Fr/Fx = (z'x/i'r)^ -1 ^ 2 - Based on 
our upper limits on the ratio, the power-law index 
of the electrons has to be < 2.7 to be consistent 
with a standard afterglow model. On the other 
hand, the case where the cooling frequency is be- 
low the optical band {y e < vr < vx), the ratio 
is expressed as Fr/Fx ~ {vx / vr) p / 2 . When we 
apply the constraint in the ratio of Fr/Fx 1$ 590, 
we find that p has to be smaller than 1.7. 

If the cooling frequency is between the optical 
and the X-ray bands [vr < v c < vx), we can 
put a constraint on the density of the circumburst 
medium (no). In this situation, the flux ratio is 
expressed as, 

F r /Fx = 14.3(n x /^)^- 1)/2 (l+ Z ) 1 / 4 £ 5 1 2 /4 €| /4 _ 2 nJ /2 , 

where t b, -2 is the fraction of energy that goes into 
the electrons in units of 10 -2 and E52 is the energy 
in units of 10 52 erg (Nysewander et al. 2009). 
Applying the flux ratio limit of Fr/Fx < 590 and 
also assuming cb » 0.01, E » fS 7l u 0 and p = 

2.1, no can be constrained as no < 0.3 cm -3 . For 
p = 2.5, we find no < 0.01 cm -3 . Therefore, we 
find a significantly smaller no compared to that of 
a typical value for long GRBs of no « 1 — 10 cm -3 
(e.g., Soderberg et al. 2006), which is likely to 
be required to explain our optical afterglow limit 
for GRB111117A. This finding is consistent with 
several other short GRBs with no optical afterglow 
detection (Berger 2010). 

We have discussed the optical afterglow limit 
based on the external shock model. However, 
there is a possibility that a significant amount 
of the optical afterglow flux was extinguished by 
the host galaxy. When we fit the X-ray afterglow 
spectrum to a power-law model with the intrin- 
sic absorption at z=1.31, the intrinsic Nh is esti- 
mated to be 7.2to'g x 10 21 cm -2 . A significance 
of the Nh detection is greater than 2 a but less 
than 3 a confidence level. Assuming a host ex- 
tinction law similar to the Milky Way galaxy, Ay 
is 4.1 mag (Predehl & Schmitt 1995). Therefore, 
a significant amount of absorption in optical flux 


is expected from the X-ray column density mea- 
surement. On the other hand, it is still not clear 
whether it is possible to have such high extinction 
at the outskirts of the host where the X-ray af- 
terglow is indicated. Moreover, the amount of ex- 
tinction which we derived from the SED fit of the 
host is Ay = 0.2 — 0.5 mag (see §6). At this stage, 
the origin of the large column density seen in the 
X-ray afterglow of GRB 111117A is still puzzling. 
Further deep optical and/or infrared imaging or 
spectroscopic observation of the host is required 
to understand the GRB environment. 

The projected offset between the afterglow lo- 
cation and the host galaxy center is 8.4 ±1.7 kpc 
using the estimated redshift of 1.31. Although this 
offset is larger than the median projected offset of 
about 5 kpc for previously studied short GRBs 
(Fong et al. 2010), it is within the offset distri- 
bution of short GRBs. Using the projected off- 
set of r = 8.4 kpc and the stellar age of r = 0.1 
Gyr, the minimum kick velocity, v = r/r, is es- 
timated to be « 80 km s -1 . The estimated kick 
velocity is similar to or possibly larger than the 
inferred kick velocity of GRB 060502B (Bloom et 
al. 2007). Using the typical age of 1-10 Gyr 
in the early-type short GRB hosts such as GRB 
050509B, GRB 070809 and GRB 090515 (Bloom 
et al. 2006; Berger 2010), the minimum kick 
velocity is estimated to be « 1-8 km s -1 . 

In this paper, we have reported the prompt 
emission, the afterglow and the host galaxy prop- 
erties of short GRB111117A. The prompt emis- 
sion observed by the Swift BAT and the Fermi 
GBM showed 1) a short duration, 2) no extended 
emission, 3) no measurable spectral lag and 4) a 
hard spectrum. All those properties can securely 
classify this burst as a short GRB. Although the 
optical afterglow has not been detected by our 
deep observations by TNG and GTC, our rapid 
Chandra ToO observation provides a sub-arcsec 
position of the afterglow in X-rays. This Chandra 
position is crucial to identify the host galaxy and 
also to measure the significant offset of 1.0" be- 
tween the host center and the afterglow location. 
Our deep near infrared to optical photometry data 
of GMG, TNG, NOT, GTC, UKIRT and CFHT 
enable us to estimate the redshift of the host to 
1.31. The observation of GRB 111117A suggests 
that X-rays are more promising than optical to 
locate short GRBs with sub-arcsecond accuracy. 
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Combining the sub-arcsecond afterglow position 
in the X-ray and the deep optical images from the 
ground telescopes, we successfully investigate the 
hcst properties of GRB 111117A deeply even with- 
out an optical afterglow. Rapid Chandra ToO ob- 
servations of short GRBs are still key to increasing 
the golden sample of short GRBs with redshift to 
pin-down their nature. 
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(2011-11-17 12:13:41.9216UT) 

Fig. 1. — The background subtracted 5 ms light 
curves of Swift BAT (15-25 keV, 25-50 keV, 50- 
100 keV and 100-350 keV) and Fermi GBM (8-200 
keY and 200-400 keV). The bottom panel shows 
the hardness ratio between the 100-350 keV and 
the 25-50 keV of the BAT data. 



Fig. 2. — Fluence ratios between the 50-100 keV 
and the 25-50 keV band versus T 90 are shown for 
GRB111117A (red) and the Swift BAT GRBs. 
The values of the Swift BAT GRBs are extracted 
from Sakamoto et al. (2011a). 



Fig. 3.— GTC r image (17" x 17") with the XRT 
90% error circle in blue and the Chandra 1 sigma 
error circle in red. 
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Table 1 

List of reference stars. 


R.A. (J2000) 

Dec. (J2000) 

R2 

[mag] 

ID of USNO-Bl 

00:50:45.254 

+23:00:59.43 

19.76 

1130-0013893 

00:50:46.923 

+23:01:47.81 

18.66 

1130-0013901 

00:50:52.923 

+23:01:19.35 

16.43 

1130-0013919 

00:50:36.987 

+23:01:25.94 

18.58 

1130-0013859 

00:50:43.442 

+22:59:06.91 

17.16 

1129-0013362 

00:50:48.207 

+23:00:36.76 

15.90 

1130-0013906 

00:50:41.781 

+22:59:33.17 

16.78 

1129-0013356 

00:50:31.795 

+22:01:36.78 

18.74 

1130-0013838 

00:50:40.927 

+23:02:53.55 

16.41 

1130-0013875 

00:50:51.577 

+23:03:10.27 

15.42 

1130-0013916 


Table 2 

Log of optical and near-infrared observations of GRB 111117A. The magnitudes are 

CORRECTED FOR GALACTIC EXTINCTION. 


Time since the trigger* 
(day) 

Telescope 

Instrument 

Filter 

Exposure 

00 

Afterglow* 

(mag) 

Host* 

(mag) 

0.083 

GMG 

YFOSC YSU 

R 

600 


> 22.1 

0.104 

GMG 

YFOSC YSU 

z 

900 


22.9 ± 0.3 (AB) 0 

0.300 

TNG 

LRS 

R 

1800 

> 24.1 

23.49 ± 0.29 

0.329 

GTC 

OSIRIS 

9 

800 


24.05 ± 0.14 (AB) 

0.342 

GTC 

OSIRIS 

r 

1200 

> 24.9 (AB) 

24.17 ± 0.10 (AB) 

0.354 

GTC 

OSIRIS 

i 

360 


23.92 ± 0.20 (AB) 

0.358 

NOT 

ALFOSC 

z 

3600 


>22.5 (AB) 

0.371 

NOT 

ALFOSC 

R 

3000 


23.20 ± 0.25 

1.5 

NOT 

ALFOSC 

R 

2400 


23.26 ± 0.22 

11.4 

TNG 

LRS 

R 

3600 


23.43 ± 0.13 

11.8 

Subaru 

IRCS 

K' 

780 


> 19.95 

14.4 

GTC 

OSIRIS 

T 

2400 


24.20 ± 0.07 (AB) 

19 - 72 f 

UKIRT 

WFCAM 

K 

8640 


20.91 ± 0.12 

42.7 

CFHT 

WIRCam 

j 

4140 


21.7 ± 0.2 


*The mean time of the stacked observations with the corresponding exposure time. 

*UKIRT data were collected at the multiple epochs between 19 and 72 days after the burst. 
*Upper limit is in 3 sigma confidence level. 

°Pos9ible host detection. 
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Flux (0.3-10 keV) Icrg cm- 2 sr 1 ] 


Table 3 

Time-integrated spectral parameters of GRB 1111 17A. 


Instrument Model 


C(GBM) 


BAT 

GBM 

GBM 

GBM 

BAT-GBM 
BAT- GBM 
BAT-GBM 


-0 eso+O- 24 
u. O^_0 22 

-1 4.4+U.Ot! 

-0 4n+ 6 - 36 

■i 07+O. U 5 
07 

-0.28in.26 

-0.28±gi| 


<-2.6 450ij2s 


420±iTg 
<-2.5 420t|J° 


0 . 62 i°;if 


49.5/57 (0.87) 
636.8/604 (1.05) 
578.9/603 (0.96) 
578.9/602 (0.96) 
729.5/662 (1.10) 
627.7/661 (0.95) 
627.8/660 (0.95) 



Time since tbc trigger [s] 


Fig. 4. — Comparison of X-ray afterglow light 
curves of long-lived and short-lived short GRBs 
observed by Swift XRT and GRB 111117A (Ted). 
The Chandra data point of GRB111117A is 
shown in red filled circle. The long-lived short 
GRBs include in this figure are GRB 050724, 
GRB 051221A, GRB 051227, GRB 060313, 

GRB 061006, GRB 061201, GRB 061210, 

GRB 070714B, GRB070724A, GRB 070809, 

GRB 071227, GRB 080123, GRB 080426, 

GRB 090426, GRB 090510, GRB 090607, 

GRB 09062 IB and GRB091109B. The short- 
lived short GRBs include in this figure are 
GRB 050509B, GRB 050813, GRB 051210, 

GRB 060502B, GRB 060801, GRB 061217, 

GRB 070429B, GRB 080503, GRB080702A, 

GRB 080905A, GRB 080919, GRB081024A and 
GRB 08 1226 A. 



Fig. 5. — Deep optical TNG (R; 1.4' x 1.2') and 
GTC (r; 1.1' x 1.0') images of two epochs. The 
right panel shows the subtracted image of the first 
and second epoch. No significant residuals are 
seen in both TNG and GTC subtracted images 
at the host location (red circle). 
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Time since the trigger [hr] 

Fig. 6. — Top: optical fluxes of the first op- 
tical detection (black circle) or an upper limit 
(filled triangle) of short GRBs (Berger 2010) is 
shown as a function of the trigger time. The 
TNG and GTC upper limits of the optical after- 
glow of GRB 111117A are shown in green and blue 
filled triangle. Bottom: optical light curves of 
GRB 1111 17A in R and r band are shown. The 
plct includes R band measurement from GMG, 
TNG and NOT, and also the r band measurement 
from GTC. 



Fig. 7. — Multi-color images at the field of 

GRB 1111 17A. From left to right, and top to bot- 
tom, the images are TNG R, NOT R, GTC g, 
GTC r, GTC i, GMG 2 , CFHT J and UKIRT K. 
The host galaxy is marked in a green circle. The 
X-ray afterglow position determined by Chandra 
is marked as a red cross. The image scale is 17" x 
17". All the images are smoothed by the Gaussian 
function with 3 pixel radius. 


Flux density (Jy) 




Fig. 8. — The SED fit to the photometric data ( g , 
r, t, 2 , J and K) using the templates of the single 
stellar populations model (Maraston 2005). The 
GTC spectrum is shown in gray. 





% 



Fig. 9. — Top: Posterior probability distribution 
of the photometric redshift by the SED fit of the 
host. All 269 SED templates axe used. Bottom: 
Posterior probability distribution of the photomet- 
ric redshift of the host by 41 SED templates of a 
solar metallicity and a red horizontal branch mor- 
phology with stellar age from 20 Myr to 15 Gyr. 



o oj i.o i.s 


Redshift 

Fig. 10. — Magnitude of the short GRB hosts 
(from Berger 2010) as a function of redshift. The 
right panel shows the magnitude of the hosts with- 
out a confirmed redshift. GRB111117A (z=1.31) 
is shown as a red star. 



log (EJ [erg] 


Fig. 11.— Comparison of E^ between short (up- 
per panel) and long (lower panel) GRBs. The 
short GRB E- xso values are from (Berger 2010) 
(samplel), and the long GRB E lao values are from 
Nava et al. (2012) (only Swift long GRBs). 
GRB111117A (3.4 x 10 51 erg) is located at the 
high end of E iao distribution of short GRBs. 
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